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We study nonequilibrium effects at the QCD phase transition within 
the framework of Polyakov loop extended chiral fluid dynamics. The quark 
degrees of freedom act as a locally equilibrated heat bath for the sigma 
field and a dynamical Polyakov loop. Their evolution is described by a 
Langevin equation with dissipation and noise. At a critical point we ob- 
serve the formation of long-range correlations after equilibration. During 
a hydrodynamical expansion nonequilibrium fluctuations are enhanced at 
the first order phase transition compared to the critical point. 

PACS numbers: 25.75.-q, 47.75. +f, 24.60.Ky, 25.75.Nq 

1. Introduction 

Presently the knowledge of the QCD phase diagram is still limited. 
While lattice QCD calculations tell us that at vanishing baryochemical po- 
tential there is a crossover [I], only model studies allow to explore the high 
baryon density regions. There one expects a first order transition at large 
fi ending at a critical point (CP) [2J. The CP is expected to be detected 
in heavy-ion collisions through event-by-event fluctuations of quantities like 
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transverse momentum or particle multiplicity OS]. Nevertheless, phenom- 
ena like critical slowing down as well as the finite system lifetime and size 
prevent the correlation length from diverging [5] which thus weakens these 
signals. It is therefore important to include nonequilibrium effects and the 
dynamics of the system to estimate the experimental signatures of the crit- 
ical end point. A promising ansatz to study the QCD phase transition in 
such a setting is provided by the framework of chiral fluid dynamics [SJ [7J [8] . 
Here, the basic idea is to propagate the order parameter of chiral symmetry 
explicitly by a Langevin equation, while the heat bath is given by a fluid 
dynamically expanding medium made out of quarks. Presently, we extend 
this model with the Polyakov loop to consider both the chiral and the de- 
confinement transition. Here, the Polyakov loop is treated as an effective 
field which is propagated by a phenomenological Langevin equation. We 
present results of temperature quenches in a box and study the evolution of 
fluctuations during the expansion of a hot plasma droplet. 



2. Chiral fluid dynamics with a Polyakov loop 

For investigation we use the Polyakov loop extended quark meson model 
[3] with the Lagrangian 

C = q [i (7^ - igqcD-fAo) -go]q+\ ( V) 2 -U{o)- U(£, I) , (1) 

where q = (u, d) is the constituent quark field, Aq the temporal component 
of the color gauge field and a the mesonic field. The pion degrees of freedom 
are neglected in the present study. The potential for the sigma field reads 

U(a) = ^{a 2 -u 2 ) 2 - h q a-U . (2) 

The temperature dependent Polyakov loop potential is chosen in a polyno- 
mial form [TO] following 

« ( t ,i) = -*f> (n> + KT) - 1 (f +?) + % {v? + KT) 2 ■ (3) 

Integrating out the quark degrees of freedom, which will constitute the heat 
bath, we obtain the grand canonical potential. At \xb = 0, I = I and in the 
mean-field approximation it is [S] 



4NfT [ ^ [l + Me- pE + ^~ WE + e" 3/3B l . (4) 

J {^r l J 



For different coupling strengths g the effective potential V e ff = U + U + Qq q 
shows different characteristic shapes at the transition temperature: for g = 
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4.7 we see two degenerate minima, see figure HJ while for g = 3.52 the 
potential is very broad and flat at the minimum, figure [2j This resembles 
the situation at a CP. These choices of the quark-meson coupling are of 
course unphysical but allow a first qualitative study of effects at several 
types of transition 0. 
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Fig. 1. Effective potential at g = 4.7, 
first order, T c = 172.9MeV. 
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Fig. 2. Effective potential at g 
3.52, CP, T c = 180.5M e y. 



In [8] the coupled dynamics for the sigma field and the quark fluid were 
derived self-consistently. The sigma field is propagated by a Langevin equa- 
tion 



d.d^a + r] a d t o + 



dV c fi 
da 



(5) 



The explicit form of the temperature dependent damping coefficient r] a 
together with the correlator for the stochastic noise field have been derived 

We allow for a dynamical evolution of the Polyakov loop by adding a 
kinetic term in the equation of motion [11] . For the full nonequilibrium 
description we also need to add a damping term r\g ~ 1/ 'fm [12] and impose 
the dissipation-fluctuation relation: 



2N n 
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d.d^iT 2 + m d t e + 



dv r 



eff 



Of 



(6) 



(UtMt')) = yS(t-t')2r k T 



(7) 



Note at this point that the Polyakov loop is originally defined only in equi- 
librium and it is not a priori clear what the correct dynamics are [10] . This 
approach is, therefore, purely phenomenological. 



Note that in principle one has to choose g such that ga reproduces the constituent 
quark mass in vacuum. This would give a value of g ~ 3.2. 
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The quark fluid is propagated via the equations of ideal relativistic fluid 
dynamic using the energy-momentum tensors of the liquid, the u-field and 
the Polyakov loop 

dr(Tf + TF' + T£ v )=Q. (8) 

3. Numerical results 

3.1. Temperature quench in a box 

We put fields and fluid in a box with periodic boundary conditions. 
The fields are initialized at a temperature above T c . Then the temperature 
is quenched to a value below the transition point and the quark bath is 
added. The fields lose energy through damping, transferring this amount 
of energy to the fluid via equation ([8]) which leads to a subsequent increase 
of the temperature, see also [13]. We find quench temperatures such that 
the system relaxes near the transition temperature. In the Polyakov loop, 
we can observe an interesting phenomenon. Figures [3] and [5] show the value 
of £ in ^-direction y = z = evolving in time for both transition scenarios 
during the relaxation process. At the CP, one observes the formation of 
long-range fluctuations over both space and time, an effect that does not 
occur at the first order transition. 




Fig. 3. Fluctuations at the first order 
transition during relaxation process, 
quench from T = 180 MeV to T = 
140 MeV. 
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Fig. 4. Long-range fluctuations at the 
CP after relaxation process, quench 
from T = 186 MeV to T = 166 MeV. 



3.2. Fluid dynamic expansion 

For this simulation we initialize an ellipsoidal region with a temperature 
of T = 200 MeV, smoothed by a Woods-Saxon-function at the edges. The 
fields and fluid are set to their equilibrium values throughout the lattice. We 
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let the system expand by full 3 + 1 dimensional fluid dynamics and measure 
the average temperature in a fixed central volume as a function of time. The 
result is shown in figure [5j While for the CP T decreases monotonically 
with time, we can observe a reheating at the first order transition as a 
consequence of the forming of a supercooled phase below the transition 
temperature that finally decays to the global minimum and transfers its 
energy into the fluid. 
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Fig. 5. Temperature of the quark fluid, for the first order scenario reheating after 
formation of a supercooled phase is observed. 




Fig. 6. Sigma fluctuations grow large 
at the first order transition compared 
to the CP. 
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Fig. 7. Polyakov loop fluctuations 
grow large at the first order transi- 
tion compared to the CP. 



In figures [6] and [7] we show the evolution of nonequilibrium fluctua- 
tions for a and I. These are defined as (Ad) = \l {{a — cr eq ) 2 ) and (A£) = 

{(£ — £ eci ) 2 ), respectively. For both order parameters we find a large in- 
crease of the nonequilibrium fluctuations in a scenario with a the first order 
transition compared to the scenario with a CP. This is again caused by the 
large deviations from equilibrium that occur during supercooling. We see in 
both figures a second smaller increase, this arises when parts of the system 
cross the transition temperature a second time after the reheating. 
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4. Conclusions 

We presented a dynamical model to study the chiral and deconfine- 
ment phase transition of QCD. The nonequilibrium evolution of both or- 
der parameters and their interaction with the quark fluid was described 
by Langevin equations. At the critical point we observed the formation of 
long-range fluctuations. During the hydro dynamical expansion of a finite 
size system our model shows nonequilibrium effects like supercooling and 
reheating of the quark heat bath at the first order phase transition. There 
we also found an enhancement of nonequilibrium fluctuations in a first order 
phase transition scenario compared to an evolution through the CP. 

This work was supported by the Hessian LOEWE initiative Helmholtz 
International Center for FAIR. 
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